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Multiplicity dependence of Bose-Einstein correlations in pp reactions: 
a discussion of possible origins* 
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The observed pronounced multiplicity dependence of correlation functions in hadron-hadron reactions and 
in particular of Bose-Einstein correlations provides information about underlying physics. We discuss in this 
contribution several interpretations, giving special attention to the string model for Bose-Einstein correlations of 
Andersson and Hofmann, as well as the core-halo picture of Csorgo, Lorstad and Zimanyi. 



1. Introduction 

It has been well known for years that 
Bose-Einstein (BE) correlations measured in 
hadron-hadron reactions at sufficiently high en- 
ergies exhibit a pronounced dependence on mul- 
tiplicity in the form of the strength parameter A. 



|4[|5[] or at least much weaker, especially with a 2- 
jet selection Heavy ion reactions, on the other 
hand, see A values decreasing with multiplicity 
at lower multiplicity densities (lower A reactions) 
but increasing againQ at higher multiplicity den- 
sities 0]. 

These findings may suggest a possible com- 
mon explanation in terms of superposition of sev- 
eral sources or strings, where each of them sym- 
metrizes separately according to the model of An- 
dersson and Hofmannl^. No BE correlations are 
then expected between decay products of differ- 



*Proceedings of the 9th International Workshop on Mul- 
tiparticle Production, Torino, June 11-18, 2000, edited by 
A. Giovannini and R. Ugoccioni; Nuclear Physics B Sup- 
plement (to be published). 

^ It should be borne in mind, though, that A measure- 
ments of different experiments are not easily compared 
because of varying systematics such as the chosen fit func- 
tion and fit range, the "cleanness" of the pion sample, 
different treatments of Coulomb corrections etc. The de- 
scribed tendencies in A are therefore more qualitative than 
quantitative in nature. 



ent strings ||]. The following arguments could 
be used: In the framework of the Dual Parton 
Model for hadron hadron and heavy ion reac- 
tions Ip^ , it is expected that higher multiplici- 
ties correspond to a higher number of strings or 
chains. This can explain the multiplicity depen- 
dence in hadron-hadron reactions. In e+e~ reac- 
tions, with a selection of two-jet events, only one 
single string is formed and consequently one does 
not expect much dependence of BE correlations 
on multiplicity. In heavy ion reactions, there is 
again an increasing number of strings with in- 
creasing multiplicity density, but eventually the 
densely-packed strings coalesce until they finally 
form a large single fireball. This picture can qual- 
itatively explain why the BE signal is decreasing 
at low A and finally increasing again [Q. A clean 
situation of a superposition of two strings occurs 
if two M^-bosons are coproduced in e+e^ reac- 
tions and both decay hadronically. If there are 
no BE correlations between the pions from dif- 
ferent strings, the A values are expected to be 
diminished. Unfortunately, these measurements 
are hampered by low statistics and experimental 
difficulties ff^-Pl- 

It should be stressed, however, that this is not 
the only possible explanation of experimental be- 
haviour. In the case of hadron-hadron reactions, 
for example, several alternative explanations do 
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exist [0, among them the core-halo picture ||l5[| 
which connects consistently the multiplicity de- 
pendence of correlation functions of like-charge 
with those of opposite-charge pion pairs. 

From all these considerations, it should be clear 
that there is important information in the ob- 
served multiplicity dependencies, in particular 
when comparing different reaction types. 

This contribution concentrates on a discussion 
of possible origins which could lead in the case of 
high energy hadron-hadron reactions to the ob- 
served multiplicity dependence. We first adress 
in Section 2 the question of the influence of jet 
production by using again pp collisions at ^/s — 
630 GeV measured in the UAl detector @,[6). 
Low-pT and high-pr subsamples are investigated. 
Using the low-pT subsample where the influence 
of jets is removed to a large extent, we finally 
discuss in Section 3 possible underlying physics. 

2. The influence of jets 

It is well known that in hadron-hadron 
reactions the probability of jet production rises 
with energy and multiplicity. The multiplicity 
dependence of correlation functions can be in- 
fluenced by this hard subprocess. To investigate 
this influence, we used a data sample similar to 
that in Q but with larger statistics: It consists 
of 2,460,000 non-single-diffractive pp reactions at 
yfs — 630 GeV measured by the UAl central de- 
tector Q . Only vertex-associated charged tracks 
with transverse momentum pT > 0.15 GeV/c and 
|?7| < 3 have been used. We restrict the azimuthal 
angle to 45° < |0| < 135° ("good azimuth"). 
These cuts define a region in momentum space 
which we call fi. 

Since the multiplicity for the entire azimuthal 
range is the physically relevant quantity, we select 
events according to their uncorrected all-azimuth 
charged-particle multiplicity N. The corrected 
multiplicity density is then estimated as twice n, 
the charged-particle multiplicity in good azimuth: 
{dNc/drf) ~ 2{dn/dr]). 

The measured quantities are K2, the second 
order normalized cumulant correlation functions 
in several multiplicity intervals N e [A, B] ; for a 
complete definition, refer to |14|. They are mea- 



sured for pairs of like-sign (is) and opposite-sign 
(os) charge separately as 



K2 iQ\AB) 



n(n-l)^. 



ri^{Q\AB) - 1,(1) 



K2' "^QIAB) ^ ^^±^ r°%Q\AB) - I , (2) 



where Wis (= ^+ — ^-) , rrfn_ and n(n— 1)^^ are 
the mean numbers of positive or negative par- 
ticles, OS pairs and is pairs respectively in the 
whole interval f2 and in the multiplicity range 
[A, B]. The prefactors in front of the normalized 
density correlation functions r correct for the bias 
introduced by fixing multiplicity |jl^. The func- 
tions Tj'* and are defined for £s and os pairs 
in the correlation integral description |jl^ 



(Q) = 



Pi®p{'{Q) 
pV{Q) 
pt®pt{Q) 

pi®pT{Q) 

pV{Q) 

pt®Pi{Q) Pi®pt{Q) 



P2~iQ) 

Pi®Pi{Q) 



P2^iQ) 



(3) 



(4) 



with Q = ■\/— (pi — ^2)^ the spacelike four- 
momentum difference of the pion pairs. 

Internal cumulants (|l]) and (^) are analysed in 
three samples as a function of pion transverse mo- 
mentum pt as follows: 

(i) An all-pT sample of all like-sign and 
opposite-sign pion pairs in 17. This is shown 
in Fig. 1 for three representative multiplic- 
ity densities. 

(ii) A low-pT subsample containing only 
charged particlesQ with p-r < 0.7 GeV/c. 
Also removed from the subsample were 
entire events containing either a jet with 
Et > 5 GeV or at least one charged par- 
ticle with Pt > 2.5 GeV/c. This is shown 
in Fig. 2 with the same three multiplicity 
selections. The prefactors entering Eqs. (|^), 



The pion sample contains about 15% kaons. 
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(H) are calculated in this case by using only 
charged particles 0.15 < pt < 0.7 GeV/c. 

(iii) A high-pT subsample where only charged 
particles with px > 0.7 GeV/c are consid- 
ered. This is shown in Fig. 3. It has been 
demonstrated previously that all high- 
Pt particles stem from jets or minijets. 

From Figs. 1-3, we see that whereas the low- 
Pt subsample behaves very similarly to the all- 
Pt sample, there is a pronounced increase in the 
strength of correlation functions for the high-p-r 
case (iii) (note the different scale on the plot!). 
This can be interpreted as the influence of jets, 
which are inherently spiky in nature. In this case, 
BE correlations are hence mixed up with corre- 
lations originating from jets, so that it would be 
difficult to measure them separately. 

We note further that £s functions in Fig. 3 re- 
veal a crossover in the region Q ~ 1 GeV, mak- 
ing the determination of the multiplicity depen- 
dence of if I highly Q-dependent. The fit to this 
high-pT sample using eq. (^) gives "radius pa- 
rameters" R which decrease with increasing mul- 
tiplicity, in contrast to the samples (i) jlj] and 
(ii). Fig. 3 shows also a pronounced secondary 
peak after a minimum at Q ~ 1.4 GeV which can 
be attributed to the onset of local px compensa- 
tion with two back-to-back particles with at least 
Pt = 0.7 GeV/c corresponding to the cut applied 
in this subsample . 

The fits performed to £s pair data in Fig.l - 
Fig. 3 are exponential, 

Ki'%Q)^a + Xe-'^'^ . (5) 

The corresponding dependence of A on multiplic- 
ity is used in subsequent figures. 

Fig. 4 compares the multiplicity dependence of 
the all-pT and low-pT^ samples in the small-Q re- 
gion (Q = 0.1 GeV). The two samples differ only 
slightly in their respective cumulants as well as 
their A values. 

Fig. 5 shows the corresponding multiplicity de- 
pendence in the high-pr sample, once again on a 
larger scale in if^- The influence of jets shows 
up dramatically: all correlation functions are in- 
creased in height, and the os functions do not 



show a multiplicity dependence for particle den- 
sities dNc/drj < 3.3. A decrease like in (i) and 
(ii) is probably compensated by increasing jet ac- 
tivity. 

The interplay of BE correlations and resonance 
production with the onset of jet production and 
the transition from soft to hard interactions are 
interesting questions in their own right and can 
be studied with samples like the high-p^ one. We 
will, however, concentrate in the following on the 
low-pT subsample. 

3. Multiplicity dependence of the low-pT 
sample 

The results of a study with the all-pT sample 
(i) have been published in Ref. |Q. In Fig. 6, we 
plot for the low-pT sample (ii) the same ratio 

X| (Q I jdNjdri ^ 6.9)) 

X| (g I {dN,/d7^ = 1.2)) ^ ' 

for like-sign and opposite-sign pairs respectively, 
while Fig. 7 shows the behaviour of the low-pr 
cumulants for fixed Q but varying multiplicity. 
These figures show that the aW-pr results found 
previously remain valid for the low-px sample 
also, namely: 

- The like-sign and opposite-sign cumulants 
have very similar multiplicity dependence 
when compared in the same Q region. 

- Cumulants behave distinctly differently at 
small and large Q: at small Q, the 
multiplicity dependence of both samples 
is weaker than 1/Nc, while at large 
Q { > 2 GeV) the cumulants are negative 
and follow roug hly a 1/Nc law|. 

- A third region around Q = 1 GeV shows 
small and rapidly changing cumulants. 

In the following, we discuss four possible explana- 
tions of these phenomena. It should be stressed 
that the LUND Monte Carlo model (PYTHIA) 
cannot reproduce the multiplicity dependence of 
correlation functions and in particular of the BE 
effect [|0). 

^ For simplicity we write Nc instead of dNc/dr] here and 
below (TVc = 6 ■ dNc/drj). 
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3.1 Bose Einstein correlations are the re- 
sult of symmetrization within individ- 
ual strings only Q : When several strings 
are produced, each string symmetrizes sep- 
arately and decay products of different 
strings would hence not contribute to BE 
correlations 

Because unnormalized cumulants of inde- 
pendent distributions combine additively, 
the independent superposition in momen- 
tum space of V equal sources/strings, each 
with a qth order cumulant Kq(pi, . . . , Pq) 
and each with some multiplicity distribu- 
tion (e.g. Poisson) results in an unnormal- 



ized combined cumulant tVq of 



(Pi 



. ,Pq) = Kg(Pl, • ■ • ,Pq) , 



(7) 



while the combined single particle spectrum 
is given in terms of individual sources' spec- 
tra by 



(8) 



so that the gth order normalized cumulants 
for V superimposed sources is given by 



V K„(pi, 



1/9 pi(pi) . . . ,Pl(Pg) 



{9) 



= A'q(pi,...,p,). 



Hence k\^^ is inversely proportional to the 
number of sources. This remains true for 
the correlation integral k'^^^ (Q) also. 

The above derivation is only for illustra- 
tion. If -fCj'^'' cx this would imply 
iiTj'^'' cx only if oc i^, i.e. for iden- 
tical sources each of fixed multiplicity. In 
reality, the assumption of equal sources is 
probably not fulfilled. The following sce- 
nario might be more realistic: "Fixing mul- 
tiplicity does not necessarily mean fixing 
the number of sources. The sources (we 
will them define below) probably do pos- 
sess a whole multiplicity distribution rather 



than a single fixed multiplicity. Our se- 
lected multipHcities range from 0.83 to 9.1, 
varying over about a factor 10. At small 
Q, however, the ifj^ vary only by at most 
a factor 3. This suggests that at the high- 
est selected multiplicity we would observe 
the superposition of only 3 sources, from 
which we are sampling their high multiplic- 
ity tails." R In the Dual Parton Model 
approach one source might be iden- 

tified with the topology of one pomeron 
exchange]^ If we select low-multiplicity 
events, we expect to select the case of 
one pomeron exchange. Multiparton col- 
lisions corresponding to multipomeron ex- 
change are expected to contribute to higher- 
multiplicity events. Estimates in ref | |2^ 
predict two- to three-pomeron exchanges 
at the highest multiplicities seen by UAl. 
The number of sources would increase cor- 
respondingly. This could explain the sup- 
pression of £s (Bose-Einstein) functions in 
Fig. 7a. However, additional assumptions 
are needed to explain the similar behaviour 
of£s and os functions in Figs. 6a and 7a and 
their 1/iV dependence at large Q ( Figs 6b, 
7b). Resonance production and colour re- 
connection effects might be candidates (see 
sect. 3.3). 

3.2 Quantum statistical approach: A 

chaoticity parameter p decreasing with mul- 
tiplicity would, in the quantum statistical 
approach 24|, decrease BE correlations 
at higher multiplicities. In this picture, 
however, the question arises what the phys- 
ical nature of the subprocess causing in- 
creased coherence at higher multiplicities 
would be. Also, the similarity of the be- 
haviour of is and os correlation functions 
in Fig. 7 is not easily explained within this 
framework. 



*This, however, implies no multiplicity dependence of cor- 
relation functions within one source, which still is a strong 
assumption. 

^ In the DPM, the exchange of one pomeron corresponds 
already to the formation of two chains or strings. We 
consider this case here as "one source". 
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3.3 BE and resonance-induced correla- 
tions combined: One could hypothesise 
that the observed multiplicity dependence 
of is correlation functions is the result of 
two processes ||25| : 

a) Bose-Einstein correlations in the classi- 
cal sense which, being a global effect, are 
independent of multiplicity [ p6| ^, and 

b) the production of higher mass-resonances 
or clusters decaying into two or more like- 
sign pions: R* tt^tt^ + X (as seen for 
example in 77' decay). If the unnormalized 
cumulants k^"'' and k^^* were wholly the 
result of resonance decays and if the num- 
ber of resonances were proportional to the 
multiplicity Nc, then oc Nc- Assum- 
ing Pi{p\Nc) oc NcPi{p) gives pi(E)pi cx 
N^, and hence after normalization, the 
resonance- inspired guess yields l/N^ be- 
haviour, 



K. 



,/ res 
I _ '^2 



Ko'"" 1 
CX 

Pi® Pi Nc 



(10) 



A mixture of processes a) and b) would give 
HQ) 



N, 



(11) 



in agreement with the behaviour of £s func- 
tions in Fig. 7a (straight hue). The 1/Nc de- 
pendence in Fig. 7b at large Q suggests the 
existence of a resonance / cluster component 
for both is and os pair production. The res- 
onance contribution to the correlation func- 
tions would be concentrated mainly at the 
region around Q ~ 1 GeV as in the case 
of p° production, which is visible as a peak 
in the os functions in Figs. 1-3. This re- 
gion is however difficult to investigate be- 
cause there the K2 are decreasing rapidly 
with increasing Q while the K2^'' are al- 
ready small. A 1 /A^c-dependence due to res- 
onances or clusters in this dominant phase 
space region around 1 GeV could presum- 
ably cause the large-Q region to follow suit 



"The observed multiplicity dependeace in ref^^pq] is com- 
pensated by our prefactors in Eqs. (hi) and (H). 



via missing pairs, thus explaining the ob- 
served dependence there. 

Once again, however, the similarity of £s 
and OS functions in the small-Q region in 
Fig. 7a can hardly be explained by assuming 
only the two components a) and b). 

One possible explanation could be the exis- 
tence of global correlations for os pairs too. 
It would give a behaviour for Ki^"" simi- 
lar to that of Eq. ( pi] ) and would explain 
its constant a{Q) part by noting that the 
number of H — pairs that can be formed 
from N-f- positive pions and N- negative pi- 
ons would be If each of these pairs 
was correlated (statistically speaking), i.e. 
if all + N-) pions were mutually corre- 
lated, the unnormalised ^2°* would be pro- 
portional to and hence after normalisa- 
tion i^l °* constant in Nc- Such a constant 
term signalizes maximum possible correla- 
tions in some events. 

More generally and following the arguments 
leading to Eq. (|l0|), we could also say: 
"If resonance production were to rise more 
quickly than cx Nc, then would de- 

crease more slowly than 1 /Nc" . 

3.4 The Core-Halo picture This pic- 

ture is currently the only one which con- 
nects the multiplicity dependence of £s 
with that of os correlation functions. The 
core-halo picture is based on the fact that 
Bose-Einstein correlations of decay prod- 
ucts from long-lived resonances are not ob- 
servable by experiments because they occur 
below experimental resolution and hence 
by definition belong to the "halo" of reso- 
nances that decay at large distances. Ex- 
amples are = 23.5 fm/c, l/F^/ = 
986.5 fm/c, l/F^ = 164400 fm/c, which 
are all unresolvable within the UAl experi- 
ment which can resolve decay products for 
distances < 6 fm only, corresponding to 
Q > 30 MeV. Because of the halo, the A 
values of BE fits are in general reduced |]15| , 



A(Q. 



fc 



-F'f,corc(p) 



(12) 
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where the second factor describes the nio- 
mentum dependence and Qmin is the smah- 
est Q value accessible by the experiment^. 
The relation between A and the fraction of 
pions emitted by the halo follows from Eq. 
( [l^ ) and from 

„ _ (^corc) _ -I (Mialo) „^ 

(iV) {N) ' 

where (A^corc) is the mean multiplicity of 
directly produced "core" pions, (A^haio) is 
the mean multiplicity of halo pions and (N) 
= (Ncorc) + (-^haio)- This meaus that 

From Eq. (jlj) it is evident that A remains 
independent of multiplicity only if (iVcorc) oc 
{N) and consequently (A^haio) oc {N). If 
however the number of halo-resonances in- 
creases faster than c>c (TV), then also the 
number of their decay particles (iVhaio) . As 
a consequence the BE parameter A will de- 
crease with multiplicity. 

A second consequence emerges immediately 
(see last sentence in sect. 3.3): The fraction 
of stemming directly from the decay 

of halo resonances will decrease less rapidly 
than (X 1/N. A previous study re- 
vealed the Q region where two-body tt'^tt^ 
decay products of the halo resonances 77, 77' 
and uj contribute, namely in 0.03 < Q < 
0.55 GeV. This is exactly the region where 
OS correlation functions indeed show a mul- 
tiplicity dependence weaker than cx as 
shown in Fig. 6b (where the 1/iV case is in- 
dicated by the dashed line). 

So far, this discussion is purely qualitative. 
How the decrease of A with Nc would com- 
pare to an effective slower-than-l/A'c de- 
crease of is, of course, a quantitative 

In is assumed that Qmin — 10 MeV and that the fit 
parameters (of Gaussian fits) are insensitive to the exact 
value of Qmin in a certain restricted region. 



question not answered by the above argu- 
ment, both because A refers to the mini- 
mum Q -value and its stability against shifts 
is not yet tested, and because the experi- 
mental fraction (A^haio) / (N) is known only 
sparsely, if at all. A more quantitative esti- 
mate has to be done in future, including the 
fact that the previously measured if^'* ^'^^ 
the whole sample are already near unity 
at Q ~ 0.03 GeV, even after subtracting the 
background contribution due to the non- 
Poissonian overall multiplicity distribution. 
Correcting for an additional halo contribu- 
tion could finally cause i^a* for Q ^ to be 
greater than 1. 

The four different explanations considered 
above each have some merit. It is clearly desir- 
able to shorten this list of candidates. We believe 
that higher order cumulants are suitable for this 
purpose: If symmetrization of individual strings 
is the right explanation, we can expect from Eq. 
(^ and Ref. |29| that the higher-order cumulants 
would decrease much faster with multiplicity than 
K2 (e.g. if I oc ). Since such a fast decrease 
is not predicted e.g. for the core- halo picture pi] ], 
the measurement of the multiplicity dependence 
of A'3 could probably decide between the two 
cases (3.1) and (3.4). 
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Figure 1. Normalised i's and os internal cumulants for the all-p^ sample (i) with multiplicity densities 
dNc/drj = 1.22 (open circles), dNc/drj — 2.72, (full circles) and dNc/drj — 6.85 (diamonds). The fits in 
a) are exponential as in eq. (^). 
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Figure 2. Same as in Fig. 1, but for the low-pT sample (ii). 
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Figure 3. Same as Fig.l, but for the high-pT sample (iii). 
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Figure 4. Comparison of the multiplicity dependence of the \ow-pT sample (ii) with that of the all-pT 
sample (i). Full lines are fits of the all-px sample to the 1/Nc behaviour as in eq. (pl|); dashed lines are 
the corresponding fits for the low-px sample. 
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Figure 5. Multiplicity dependence of the high-pT sample (iii). The lines (full lines for £s pairs, dashed 
line for os pairs and for l/{dNc/dr]) > 0.2), pro-forma fits using (Eq. (pi])), are clearly an inadequate rep- 
resentation of the data and are hence intended only to provide comparison to 1/Nc behaviour. Note that 
these high-pT cumulants exceed one by a considerable amount; high-pT data is hence clearly dominated 
by processes other than BE correlations. 
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Figure 6. The ratio of two if^ (Q) corresponding to two selections of dNc/drj as indicated, for the low-pr 
sample (ii) . The dashed line indicates the value of the ratio for the case that K^{Q) oc 1/Nc 
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Figure 7. a) Multiplicity dependence of Ki^^'' (filled circles) and (open circles), both at Q = 0.1 

GeV, b) as in a) but for the large Q = 7 GeV-region. Note from Fig. 1 that the cumulants are negative 
at large Q. For better comparison of the respective dependencies on dNc/drj, the absolute values have 
been scaled by constant factors. Dashed lines are best fits using Eq. (|Tl|). 



